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1 Introduction. 

O 

I/"") ■ Just at the beginning of quantum stochastic calculus (QSC), Hudson and Parthasarathy proposed a 

quantum stochastic Schrodinger equation linked to dilations of quantum dynamical semigroups |fl ^[ [ll]] . 
Such an equation has found applications in physics, mainly in quantum optics, but not in its full generality 
H [jj |^| . It has been used to give, at least approximately, the dynamics of photoemissive sources such 
as an atom absorbing and emitting light or matter in an optical cavity, which exchanges light with the 
surrounding free space. But in these cases the possibility of introducing the gauge (or number) process 
in the dynamical equation has not been considered. In this paper we want to show, in the case of the 
simplest photoemissive source, namely a two-level atom stimulated by a laser, how the full Hudson- 
Parthasarathy equation allows to describe in a consistent way not only absorption and emission, but also 
the scattering of the light by the atom. 

Let us recall the Hudson-Parthasarathy equation; this is just to fix our notations, while for the proper 
mathematical definitions and the rules of QSC we refer to the book by Parthasarathy fill . We denote 
by T :— T{X) the Boson Fock space over the Hilbert space X := Z ® L 2 (R+) ~ L 2 (R + ; Z), where Z is 
another separable complex Hilbert space. Let {e,, i > 1} be a c.o.n.s. in Z and let us denote by A,(t), 
A\(t), Aij(t) the annihilation, creation and gauge processes associated with such a c.o.n.s. We denote by 
E(h), h € X, the exponential vectors in T with normalization ||i?(/i)|| 2 = exp{||/i|| 2 }; E(0) is the Fock 
vacuum. We shall also use the Boson Fock spaces Tt '■= T(L 2 ([0, t]; Z)j and T l := T(L 2 ((t, oo); Z)j , for 
which we have T — Tt® T 1 , and the Weyl operators 

WtC/J^expj^jr^/jWdAtW-^cLljWjJ, /£L 2 JK + ;2). (1) 

Let H be a separable complex Hilbert space (the system space) and let H , i > 1}, {<SV/, i,j > 
1}, be bounded operators in H such that = H n , * s strongly convergent to a bounded 



operator, and J2i j &ij ® l e «)( e jl =■ S e U(H •Si Z) (unitary operators in H ® Z); we set also K, 



o 



Hq — ^ Y^j Then ( |il]| Theor. 27.8 p. 228) there exists a unique unitary operator-valued adapted 

process U(t) satisfying U(0) — 11 and 



dU{t) = {Y, R0 3 dA l ( f ) + - *«) dA ^ ^ ~ J2 R i* S v dA J (*) ~ iK ° dt ) 

V A A A A A J 



U(t). (2) 
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Now, T is interpreted as the Hilbert space of the electromagnetic field; A^{i) creates a photon of type 
j in the time interval [0, t], Aj(t) annihilates it, Ajj(t) is the selfadjoint operator representing the number 
of photons of type j up to time t and N(t) := . Ajj(t) is the observable "total number of photons up 
to time t" . We shall see in Section 3 how to choose the one-particle space Z [eq. (0)]. 

In order to describe a two-level atom, we take Ti = C 2 ; then, to fix the model, we have to determine 
the atomic operators Ho, R®, Sij on the basis of physical considerations. Let us note that, if R® — 0, the 
flux of incoming photons turns out to be equal to that of outgoing ones, but in general, for an arbitrary 
choice of the system operators, the flux conservation does not hold. Physically, flux conservation is to 
be expected when the possible processes are absorption/emission and elastic scattering; note that in the 
presence of absorption/emission this conservation cannot be instantaneous, but only in the average over 
long times. In the next section we shall require a weak form of photon flux conservation, namely only 
in the mean and for large times [eqs. (||), fll4|)]. This suffices to determine the structure of the atomic 
operators [eqs. (fTo|), (EBJ)] and to eliminate anelastic scattering; an interesting balance equation [eq. (|25|)] 
is obtained as a byproduct. We end Section 2 by studying the large-time behaviour of the reduced atomic 
density matrix [eqs. Q)-([I||)]. 

In Section 3, in order to give an example of physical consequences of the model obtained in Section 
2, we study the differential and total cross sections for the scattering of laser light by the atom, as 
a function of the frequency of the stimulating laser [eqs. (|6l]), (|62])]. The resulting line-shape is very 
interesting. Not only a Lorentzian shape is permitted, but the full variety of Fano profiles pp. 61-63) 
can be obtained [eq. (p6[)]. Moreover, the dependence of the line shape on the intensity of the stimulating 
laser is computed; in particular, the resonance position turns out to be intensity dependent [eq. ([671) 1, a 
phenomenon known as "lamp shift" |J. 

2 Master equation and flux conservation. 

First of all we want a model for an atom stimulated by a laser; this means to choose as initial state 
G H <S> T a generic state for the atom and a coherent vector for the field [Q], i.e. 

vf = C®p|^, ^H, ||£||=1, /eL 2 (K+;Z). (3) 

Moreover, we shall consider only adapted observables X t £ £(Tt ® Ft) for which we have 

(U(t)*\X f U(t)V) = (U(t)£ ® E(0)\W t (f)*X t W t (f) U(t)t ® E(0)) , (4) 

U(t):=W t (f)*U(t)W t (f). (5) 

This implies that in all physical expressions we can take / S L 2 c (R+; Z) and indeed, in order to describe 
monochromatic coherent light, we choose 

f(t) =e~ iuJt A, XeZ, wel. (6) 

Let us recall that the atomic reduced statistical operator g\(t) is defined by the partial trace 

g x (t) := Tr V {U(t)\*)(*\U(tn = Ti V {ft(f)|£ ® E(0))(t E(0)\U(t)*} . (7) 

Moreover, the quantity 

(N(t)} := (U(t)V\N(t) U(t)9) = (U(m ® E(0)\W t (f)*N(t)W t (f)U(t)( E(0)) (8) 

represents the mean number of photons up to time t, after the interaction with the atom, while (^\N(t) 
is the same quantity before such an interaction [jjj . By the theory of quantum continuous measurements 
H |l| the probability law of the counting stochastic process associated with the observables N(t), t > 0, 
could be obtained; however, in this paper we shall need only mean values such as (g|) and not the full 
theory of continuous measurements. 
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In order to formulate physical requirements, let us start by considering the case when no photon is 
injected into the system, i.e. A = 0. In these conditions it is natural to ask that the atom can emit or one 
or zero photons depending on the atomic initial state; moreover, we ask the final state to be independent 
from the initial one; this is done in the next proposition. 

Proposition 1 We assume that 

A = V£,Vt, (JV(t))<l, eaW^Ql; (9) 

moreover, we take as canonical basis {|+), |— )} in H the basis which diagonalises (r^, so that we can 
write gr = pP+ + (1 — p)P- for some p in [0, 1]. Then, apart from an exchage of roles between the two 
states |+) . |— ), we obtain g^ = P_ and 

H = -w a z , w 6M, F& = (e j \a)<T- i «e2,a/0. (10) 



Let us recall that a z , a± are the Pauli matrices, which are given by 
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(11) 

moreover, the two orthogonal projections P± can be written as P + = i(ll+cr z ) = er + cr_, P_ = |(I1— <J Z ) = 

<7-<J + . 

Proof. By using the rules of QSC one obtains that go(t) satisfies the master equation 

^p.=£ Q [e Q (t)] > £o[Q]:=-i[H ,Q} + ~Y l ([ R h>^*] + [ R he R T})> ( 12 ) 

3 

and that 

(N{t)) = j^Tt n \Y^Rf R° jgo (s)}ds. (13) 

From cqs. ( |l3| ) and (^), we obtain the necessary condition £\ Tr« {-^i*^?^,} = By *he cyclic 
property of the trace and the positivity of gP and R ( jg® Rf* , we get that this condition is equivalent to 
^f?° q = 0,Vj. 

Now, let us set i?° = Xji + yjC z + Zj<j+ + Oj<7_ (every operator on C 2 can be written in this way). 
Then, RjQ^ = gives p(xj + yj) = 0, (1 - p)(xj - yj) = 0, (1 - p)zj = 0, paj = 0. For p E (0, 1) 
this system gives R® — 0, which is not acceptable because there is not a unique equilibrium state. For 
p = we get Xj = yj and Zj = 0; we need also l^jl 2 7^ to have decay to an equilibrium state. We 
do not consider the case p = 1, because it is analogous to the previous one, apart from the exchange of 
|+) and |— ). Therefore we have R® — a.jO- + PjP+- From eqs. @ and (|l^) we have also £o[f?°J = 0, 
which reduces to [Ho, g® ] — 0. Because Hq is selfadjoint and defined up to a constant, we can write 
Hq = ia>oc z , luq £ M. Finally, by solving the master equation and computing (N(t)) from eq. (|l3j), we 
get that (iV(t)) < 1 implies = 0. ■ 

Now, let us consider the case A ^ 0. If the possible physical processes are absorption/emission of 
single photons and elastic scattering, we expect that for large times the mean number of injected photons 
(ty\N(t) = \\X\\t be equal to the mean number of outgoing photons, so we require the flux conservation 
in the mean, as done in the next proposition. 



Proposition 2 We assume eq. and 

lim /XI?,? T v = 1 . VA e Z , Vw . (14) 

Then, we have 

S = P+®S + +P-&S- , 5 ± eW(Z). (15) 
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Proof. This proof relies on some cumbersome computations; we only try to sketch it. First of all by 
QSC we obtain 

d U(t) = | £ J#(t) <^J(*) + E - *«) dA «(*) - E <^i(*) - dij C/(t) , (16) 

i#(t) := ( ej \a + F-f(t))a_ + (e 3 \F+ f(t))a+ + (e,\(S+ - l)f(t))P + + (e^S" - l)/(t))P_ , (17) 

K\{t) := Ho - i||A|| 2 + i ((A|S+A) - (al^/W) - i||a|| 2 ^ P+ + i(A|S~A) P- + 

+ i({\\F+\) + {F-\\\)-{a\S-f(t)))a++i{f(t)\a)a-. (18) 

We have used the notation 

S =: P+ ® S+ + P- (g> S~ + a + (g> F+ + a- (g> F- , S ± , F ± e C(Z) ; (19) 

the unitarity of S implies some simple relations among , F^ . 
Then, we introduce the quantity 

Y(t) := (N(t)) + \ Tr„ {a zQx (t)} - ± (^[^0 - £ ll/WII 2 ds , (20) 

whose derivative can be computed by QSC from eqs. (@), (|), By using the unitarity of S 1 and 

the fact that Y(0) = 0, we can write 

Y(t)= [\- e X+(s)\\F-\\\2 + 2Rc e X-( S ){F+\\S + \) + g - x -( S )\\F+\r}d S , (21) 
Jo 

where g\ + (t) , . . . are the matrix elements of Q\(t). 



Again by applying the rules of QSC to eqs. (uj\ and (|16|), we obtain the master equation 
dg x (t) 



dt 



-iK x (t) Qx (t) +ig x (t)K x (t)* + ^R$(t) Qx (t)R>>(ty ; (22) 



now, the Liouvillian is time dependent. In general an equilibrium state does not exist, but, for large t, g x (t) 
can contain only constant and oscillating terms; therefore, the Cesaro limit lim^oo j J Q * g x (s) ds =: g 
exists. By dividing eq. (J2l|) by ||A|| 2 i and taking the limit, we obtain 



lim , T !!l!!L - 1 = |[F+A|| 2 + 2Re^-(F+Al5+A) - g^ + (\\F + xf + A|| 2 V (23) 

where A := A/||A||. We want the r.h.s. of eq. ( p3| ) to be zero for every A, but we know that g — > P_ for 
||A|| — > and, in particular, we have g+ + — > and g A — > 0. Therefore, the flux conservation implies 
||F+A|| = for every A, i.e. F+ = 0. 

If Z were finite dimensional, the unitarity of S and F + = would imply directly F~ — 0. In the 
infinite dimensional case, one can only conclude that F~*F~ =: Q is an orthogonal projection. Some 
more considerations are needed to have F~ = 0. The line of the proof is the following. For every A with 
a non vanishing component in the range of Q the flux conservation would imply g ++ = and so g = P_ . 
Moreover, P_ would be an invariant state for the master equation (|2^); but to impose this gives F~ = 0, 
contrarily to our hypothesis. So eq. ( |l5| ) is proved. ■ 

From now on we assume eqs. (||), (|^), jTo|), ( |l5| ) to hold and, always for physical reasons, we take 

u>a > , lj > . (24) 
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Let us stress that as a byproduct of the asymptotic flux conservation (jlj) we have obtained the balance 
equation 



(N(t)) + ±Tr n {a z (g x (t) - g x (0))} = \\X\\ 2 t, 



(25) 



which means that the mean number of outgoing photons plus the mean number of absorbed photons (the 
second term reduces to g\ + (t) — g\ + (0)) is equal to the mean number of ingoing photons. 
Let us end this section by studying the asymptotic behaviour of g\(t). By setting 



g x (t) := exp{kja z t/2} g x (t) exp {-iu>a z t/2} , 
we obtain, by QSC, the master equation 



(26) 



with 



-g x (t) =£ A [g x (t)}, 



C\[q] ■= -i[H x , g] 



Rj , P-Rj * 



(27) 



(28) 



:= (e,|a)a_ + (e,|(S+ - 1)X)P+ + (e^S- - H)A)P_ , 



(29) 



H x := - (ojq - u)a z - Im (X\S + X)P+ - Im (A|S'"A>P_ 
+ ~((S- + l)\\a)a--~(a\(S- + l)\)<r+- 



(30) 



The general master equation for a two-level system is studied in [ fTo| ; in the following we shall use 
similar techniques, apart from a different parametrization of the statistical operator which turns out to 
be more convenient in our case. By setting 



Qx(t) 



u(t) v(t) 
v(t) 1 - u(t) 



< u(t) < 1 

u(t) >u 2 {t) + \v(t)\ 2 



(31) 



where the conditions on the right express the fact that g x (t) is a statistical operator, we obtain from the 
master equation 



where 



— u(t) = Gu(i)-w, 
at 



u(t) := (u{t),v{t),v(t) 



(0,(a|S-A),(S-A|a)) T , 



(32) 



(33) 



/ -Il«ll 2 
(a\(S+ + S-)X) 

V ((S+ + S-) X\a) 



-(S-\\a) -(a\S-\) \ 
-b 
-b J 



(34) 



b := - x 2 - i (Alo - Im {S + X\P a S- A)) , 



(35) 
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x 2 := ^ + \(a\AS\)\ 2 , fi 2 := ||a|| 2 + ||P ± AS*A|| 2 , AS := S + - S~ , (36) 

Alu := lu - (uj a + lm{S + \\P ± S-\)) , (37) 

a:=-^--, P a := \a) (a\ , P ± :=l-P a . (38) 
\a\ 



Moreover, we have 



detG = -||a|| 2 (Aw) 2 +r 2 /4 , (39) 



with 

T 2 := >f 4 + 4^ 2 Re(5-A|P Q (S* + + 5-)A)-4(lm(5 + A|P Q S-A)) 2 = (40) 
= (m 2 + 2|(5|S'"A)| 2 -2Re(S* + A|P Q 5"A)) 2 + |(a| (S+ + S~) A)| 2 (2fi 2 + |(5|AS*A)| 2 ) . 

Let us note that ||a|| > implies det G < and T 2 > 0. Finally, the equilibrium state is given by 

u(oo) v(oo) 

lim g x (t) { I , (41) 



t — > + oo 



v(oo) 1 — u(oo) 



where u(oo) and v(oo) are computed by equating to zero the time derivative in eq. (|32|); then, we have 
u(oo) = G _1 w, which gives 

x 2 \\P a S~\\\ 2 

u{oo) = (AWTfyl ' (42) 
v (°°) = ~ (A iy+rV4 (j + iAL, + ihn{s+ ^ PaS ~ x) ) ■ (43) 

Quantities like cjq, <5' ± are phenomenological parameters, or, better, they have to be computed from 
some more fundamental theory, such as some approximation to quantum electrodynamics. The whole 
model is meaningful only for uj not too "far" from luq and ujq must include the Lamb shifts. In the final 
results one can admit a slight w-dependence in the elastic scattering matrices S . 

3 Cross section. 

The approximations which allows to describe the electromagnetic field as a Boson field in our Fock 
space T are known as quasimonochromatic paraxial approximation (see [Q and references therein). In 
particular in this approximation the fields behave as monodimensional waves, so that a change of position 
is equivalent to a change of time and viceversa. Moreover, we do not take into account the polarization 
degrees of freedom. Therefore, the space Z has to contain only the degrees of freedom linked to the 
direction of propagation, so that we can take 

Z = L 2 (T,d 2 cr) , T = {0 < $ < 7T, < ip < 2vr}, d 2 a = sintf dtf dip . (44) 

In this section we want to compute the cross section, when the stimulating laser is well collimatcd 
and the intensity of the light is detected in directions different from the direction of propagation of the 
laser beam. So, to have a laser beam propagating along the direction $ = 0, we take 

X = V e iS \, V >0, *G [0,270, W<p)= Aq * [0 '** lW =f ; (45) 

Aw- v /27r(l — cos Atf) 



in all the physical quantities the limit Ai9 J, will be taken. Note that ||A|| = rj/A-d, because we need a 
not vanishing atom-field interaction in the limit. 

Let us consider a spherically symmetric atom. Then we have 

S=l/V2ff, S ± a = e [s± a, S ± *a = e- is± a, (46) 

where s+ and s_ are the s-wave phase shifts for the elastic scattering in the up and down atomic states 
respectively. Moreover, we set 

g ± (d):=^ Q (p ± (S ± -l)x)(^ i <p), Ag := g+ - g_ , As := s + - s_ ; (47) 

the fact that g+ and g_ do not depend on (p is again due to spherical symmetry. Now, in the limit Ai9 J, 0, 
we can compute the various quantities introduced in the previous section; by using (5 | A) = 1/2 and 
(a | {S ± — 11) A) = (e is ± - 1) /2, we obtain 

" 2 =M 2 + ^ 2 [l-cos(A S )], /, 2 = ||a|| 2 + ?? 2 ||A.g|| 2 , (48) 

r 2 = // + 2 M V + 1 r? 4 [1 - cos( As)] , (49) 
Auj = lu - [uj + if Im (g+\g-) + \AV Im Ag(0)\ , (50) 

X 2 T] 2 /4 

M(oo) ^ (A^ + rV4 ' (51) 



v(oo) — — — e 



1 „i(5+ s _) \\ a \\V 



(Aw) 2 + r 2 /4 I 2 



2 

Acj sin(As) 

4 v y 



(52) 



By direct detection, it is possible to measure the intensity of the light (or to count the photons) 
propagating in a small solid angle AT around some direction (??', </?'), different from the direction $ = 0. 
The observable "number of photons in AT up to time i" is represented by 

X t ^^(e.lW^AyW, ( ei \lsT e,) ~ e. ( (tf', <p>) e 3 {$' , iff) sin^d^'d^. (53) 

Because of Iat A = 0, which expresses the fact that the direction of detection is different from the beam 
direction, we have Wt{f)* XtWt(f) — X t . By introducing the operator "number of photons up to time t 
per unit of solid angle" 

n(0,<p;t) :=^e i {§, V >)K ij {t)e j {§,^), (54) 
i,3 

by the previous remark and eq. (Q), its mean value for •& ^ is given by 

(n(0, tp; t)) = (U t £ O S(0) | n(tf, i)*7 t £ ® £(G)) . (55) 
Then, a natural definition of differential cross section is 

where A is a kinematical factor to be determined and with dimensions of an area. To determine A let 
us consider the cross section for purely elastic scattering, for which the Bohr-Peierls-Placzek formula (or 
optical theorem) gives cr el ($, ip) = |g($)| 2 , <J^ OT — 2 ^ Im g(0); the total cross section is the integral of the 
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differential one on the whole solid angle. In our case we have <7 cl ($, ip) = Aq |((S f± — 11) A) <p)\ 2 /|| A|| 2 
and, by the unitarity of S ± , 



|A|| 2 



| ^±-1)^1*= - 



lmie~ iS ((S ,± - 11) A) (0). 



Then, by imposing the optical theorem, we get Aq = ( 2 fp) 2 ^aq)' 1 ' ^° now we nave n °t taken into 
account the polarization degrees of freedom. If they are taken into account and the cross section for not 
polarized light is considered, a 3/2 extra-factor is obtained (0 pp. 532-533) and eq. (Eft) becomes 



3 1 

lim - (n($, (f; t)) 



By differentiating eq. (p5J), we obtain 



(57) 



- (n(tf, tp; t)) = Tr„ {Rtf, <pf R{d, <p) g x (t)} , 



(58) 



R{$, tp) := a(0, <^)a_ + ((S+ - 11) A) (tf, + ((5" - 1)A) (0, ^)P_ ; 
moreover, from eqs. ( |58| ) and (|57|) we get 

_ /2vrc\ 2 3 



ct(0 )¥ >) 



\ u> J 2irri 2 



(59) 



(60) 



Finally, by eqs. (|4l|), (pl[), (p2|), (|60[), we obtain the differential cross section and, by integrating it, the 
total one: 



2ttc 



27T77 2 \ U) J { \ 47T 

1 



1 



V 



9-W 



40F 



1) 



r?||a|| 



Ree~ iS v(oo) 



9-W 



9 + m + ^ (e- - 1) 



(l - u(oo)) + 

1 / 



u(oo) + 



40F 



1 



(61) 



3 /2ttc 
2vr l~uT 



| 2 + J (l-coss_) 



(l — u(oo)) + it(oo) X 



(62) 



1.9+1 



1 . 

- (1 -coss^ 



COS S- — 



2||a|| 2 Aw 



sin s_ 



2^ 2 



sin(As) sin s_ 



Let us end with some comments about er TOT . When the elastic scattering is negligible, i.e. when 
g± = and s± = 0, eq. (|62]) reduces to 



_3_ /2ttc\ 2 T 2 /4 



1 



2tt V w / (Aw) 2 +r 2 /4 l + 2r; 2 /||a|| 2 ' 



(63) 



with r 2 = ||aj| 4 + 2||a:|| 2 7/ 2 and Alu = ui — ojo. For a laser with negligible intensity, i.e. when 77 J, 0, 
eq. ( |63| ) reduces to the cross section for resonant scattering, given in Q pp. 530-533; for r\ 7^ 0, we have 
a power broadening of the resonance line, which maintains a Loretzian shape ^). 
To simplify the general case, we set 



2Aa 



.4 



15-1 



(1 — cos s_) , 



c 



sin s_ 



(64) 
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B := 



9 2 

r 2 



r 2 



1 i 

— cos s_ cos s 4 

2 2 4 



>£ 2 cos S- + — sin( As) sin s_ 



Then, we can write 



27T / u \ 



2 

3 V2ttJ CTt 



= A- 



B + Cx 
x 2 + 1 



(65) 



(66) 



the positivity of c TOT is equivalent to A > and A(A + £>) > C 2 /4 or A = 0, B > 0, C = 0. According 
to the values of the various coefficients and mainly to the signs of B and C, different line shapes appear, 
which are known as Fano profiles (Q pp. 61-63). These shapes are typical of the interference among 
various channels, when one of them has an amplitude with a pole in the complex energy plane; in our case 
the channels are elastic scattering in the up state, elastic scattering in the down state and fluorescence. 
Whichever be the line shape, there is a strong variation of the cross section for oj around ojq + e, 



e := rj 2 Im (g+\g-) + V^V 2 Im Ag(0) 



(67) 



The intensity dependent shift s of the resonance frequency has received various names in the literature; 
a very suggestive one is lamp shift, a name suggested by A. Kastler in Note that in our two-level 
system the lamp shift is not vanishing only if the two states respond differently to elastic scattering; 
moreover, only the not s-wave contribution does matter. Let us stress that also the line width T and the 
whole line shape are intensity dependent. 
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